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We present preliminary results for Jb and masses of low- lying heavy-light mesons in the static limit. Cal- 
culations were performed in the quenched approximation using multistate smearing functions generated from a 
Hamiltonian for a spinless relativistic quark. The 2S-1S and IP— IS mass splittings are measured. Using the 
IP-IS charmonium splitting to set the overall scale, the ground state decay constant /s, is 319 ± ll(stat) MeV. 



1. INTRODUCTION 

Lattice gauge calculations of heavy-light me- 
son structure are of both theoretical and phe- 
nomenological interest. yj] One immediate goal of 
these calculations is to obtain precise quantita- 
tive results for masses, decay constants, and form 
factors in the static approximation, where the 
heavy quark propagator is replaced by a time- 
like Wilson line. One difficulty which plagued 
early, exploratory calculations of the pseudoscalar 
decay constant fg was the problem of isolating 
the ground state contribution to the propagator 
of the local weak current. Because of the prox- 
imity of excited states and their sizeable over- 
lap with the local current, a large separation in 
time was required, with an accompanying loss 
of statistics. Recent attempts to overcome this 
problem have employed nonlocal Qq operators (in 
a fixed gauge smeared over a cube|[| or wall 
source By measuring the asymptotic behavior 
of both the smeared-smeared (SS) and smeared- 
local (SL) propagators, one can reduce the sys- 
tematic error associated with excited state contri- 
butions. However, it is likely that such smearing 
functions are too crude to obtain accurate val- 



ues of the parameters of the low-lying heavy-light 
states. H This is illustrated in Table 1. 



Table 1 

The overlap amplitudes between various cube 
smearing functions and the approximate wave- 
functions for the ground state |0) and first three 
excited states (|1), |2), and |3)) for j3 = 5.9 and 
k = 0.158 on a 16 3 lattice. 



cube 



|0> 



|1> 



|2> 



|3> 



point 
3 
5 
7 
9 

11 

wall 



0.208 
0.568 
0.764 
0.800 
0.746 
0.663 
0.459 



0.234 
0.480 
0.304 
-0.024 
-0.330 
-0.546 
-0.687 



0.285 
0.362 
0.010 
-0.296 
-0.301 
-0.099 
0.477 



0.344 
0.126 
-0.220 
-0.197 
-0.122 
-0.207 
0.002 



'Based on talks presented by E. Eichten and B. Hill 



It is important, therefore, to develop new tech- 
niques which allow the extraction of the prop- 
erties of heavy-light states from relatively short 
times. Here we report preliminary results ob- 
tained using the multistate smearing method dis- 
cussed elsewhere |6|Q]. The hallmark of a pure, 
isolated ground state meson is an effective mass 
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plot which is constant in time. In Figs. 1 and 2, 
we show our results for both the SS and SL local 
effective mass plots at /3 = 6.1, k — 0.151, on a 
set of (48) 24 3 x 48 lattices. On the horizontal 
axis is the time T in lattice units. On the vertical 
axis is ]n{Css(T-l)/Css(T)). 

The SS local effective mass reaches its asymp- 
totic value around T = 2, while the SL propaga- 
tor is nearly asymptotic after T = 3. 

The results exhibit a single consistent plateau 
at ma — 0.619 ± 0.006 over a large range of T for 
both SS and SL propagators. The determination 
of the fitted value of the effective mass will be 
discussed in Section 4. We note here that the 
effect of the off-diagonal entries in the covariance 
matrix of the data for the two correlators has been 
included. @ 

These mass plots convincingly demonstrate the 
effectiveness of our smearing method in isolating 
the ground state. 
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Figure 1. k = 0.151 SS effective mass plot. 
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Figure 2. k = 0.151 SL effective mass plot. 



In the analysis presented here a two state 
smearing matrix was used. After a few iterations 
of the multistate smearing method a rough value 
of the mass parameter fi (in the SRQM) was de- 
termined. The output approximate ground state 
wavefunction at each k value |0) was extracted 
and used with the first excited state |1) gener- 
ated from the SRQM. These were the two states 
used as smearing functions 

Two points need to be made about this varia- 
tion of the multistate analysis: 

(1) The states are not exactly orthogonal. If 
x = (0|1), then an orthonormal pair of states can 
be obtained by replacing |1) by |1'} given by: 

ii') = (ii>-z|o>)A/r^ (i) 

This effect is included here but was not included 
in previously presented results. The effect on fg 
is approximately 12 percent. 

(2) The state |0) varies with light quark hop- 
ping parameter k and this variation has some sta- 
tistical fluctuation, which is reflected in the vari- 
ation with k of our results. 



2. WAVEFUNCTIONS 

The basics of the multistate smearing method 
and the choice of the Hamiltonian to gener- 
ate the smearing wavefunctions are reported by 
Thacker [Q . The agreement between the spinless 
relativistic quark model (SRQM) wavefunctions 
and the lattice QCD wavefunctions is discussed 
in detail there. 



3. RENORMALIZATION 

In order to extract continuum results from our 
lattice calculation, the matrix element calculated 
in the lattice effective theory must be related 
to the corresponding quantity in the continuum. 
This matching is done in perturbation theory in 
two steps. The first step is to relate the operator 
(in this case the axial current) in the continuum 
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effective theory to its counterpart in the full the- 
ory @@: 

We use the 1992 particle data book average jy] 
for as = 3 2 /47r at the scale /i = a -1 and mass 
rrifc = 5 GeV. In the second step the lattice cur- 
rent is compared to the corresponding current in 
the continuum effective theory P|lC|]: 

Zut = 1 + 7^20.37. (3) 

The renormalization of the axial current is then 
Z^ 1 = Z c ffZ\ at - Table |^ lists the renormalization 
constants for all the lattices analyzed here. We 
use the results for the IP- IS splitting in charmo- 
nium to obtain a -1 

Table 2 

The current renormalizations for different lattices 



lattice 





Z c s 


Z\ a t 


Z A 


12 3 x 24 


5.7 


0.90 


1.33 


0.83 


16 3 x 32 


5.9 


0.96 


1.31 


0.79 


24 3 x 48 


6.1 


1.00 


1.30 


0.77 



4. ANALYSIS 

We now discuss the statistical analysis leading 
to our results for Jb- In the first and second sub- 
sections below, we will explain the methodology 
leading to the central values and statistical errors 
quoted in our abstract and conclusions. In the 
third subsection, we will look at variations in the 
analysis procedure and their impact on these re- 
sults. Variations are partly attributable to the 
remaining systematic errors in the computation, 
allowing an estimate of the magnitude of these 
errors. 

We present results from three different lattices 
as listed in table ||. The gauge configurations 
are separated by 500 (12 3 x 24, = 5.7), 2000 
(16 3 x 32,(3 = 5.9), and 4000 (24 3 x 48,/? = 6.1) 
pseudo-heat-bath sweeps respectively. They were 



fixed to machine accuracy in Coulomb gauge us- 
ing Fourier acceleration. We use the Wilson ac- 
tion for the light quarks. We present results 
from a single k value at = 5.7 and = 6.1. 
At = 5.9, the range k = 0.154 - 0.159 in- 
cludes 0.154, 0.156, 0.157, 0.158, and 0.159. 

Table 3 

The lattices and parameters 



lattice 





confs. 


flight 


12 3 x 24 


5.7 


48 


0.165 


16 3 x 32 


5.9 


48 


0.154-0.159 


24 3 x 48 


6.1 


48 


0.151 



4.1. Effective Masses and Selection of Fit 
Range 

The effective mass plots for the 24 3 x 48 lattices 
shown in Figs. 1 and 2 exhibit the asymptotic 
ground state signal starting at quite small times. 
However, at small times, even with carefully cho- 
sen wave functions, there are significant contribu- 
tions to smeared source-smeared sink (SS) and 
smeared source-local sink (SL) correlators from 
higher energy states with the same quantum num- 
bers as the B meson. This is apparent below time 
3 in the = 5.9 n = 0.158 SS effective mass plot, 
depicted in Figure ||[ 
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Figure 3. k = 0.158 SS effective mass plot. 



From time slice 3 on, the local effective mass 
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is consistent with a constant function of time. 
This is plausible from inspection of the plot itself, 
which has the rms (not jackknife) errors plotted 
in addition to the means. To confirm this, we 
fit the logarithm of the SS correlator to a linear 
function over the range of times 2-9, and found 
X 2 =8.9 for 6 degrees of freedom. If the lower 
limit of the fit range is reduced to 1, x 2 increases 
to 12.1 for 7 degrees of freedom. These values of 
X 2 are obtained from a single fit using 48 gauge 
field configurations. To obtain the effective mass 
and its statistical error, this fit is jackknifed us- 
ing 12 subensembles, each with 4 of the 48 lat- 
tices removed. Using the fit range 2-9, we find 
ma=0.664 ± 0.011 (statistical). 

We now turn to the corresponding SL effective 
mass plot, depicted in Figure |L This correlator 
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Figure 4. k = 0.158 SL effective mass plot. 



has much smaller fluctuations (note the change in 
scale for the ordinate), and is capable of reveal- 
ing statistically significant variations in the local 
effective mass over the range 2-9. Such varia- 
tions are apparent in effective mass plots at lower 
values of k (not depicted), and can be ascribed 
to an admixture of higher energy states with the 
same quantum numbers as the B meson. Calcu- 
lations of x 2 f° r various fits at k — 0.154 suggest 
that the lower limit of the fit range should be 
increased to 5. Using the fit range 5-9, for the 
k=0.158 SL correlator we find that ma = 0.667 
±0.009 (jackknife). 



The n — 0.158 effective masses from the two 
correlators appear to be consistent to within one 
standard deviation. However, monitoring x 2 for a 
simultaneous fit that requires the two correlators 
to have a common effective mass gives a more pre- 
cise indication, since the statistical fluctuations of 
the SS and SL correlators may be correlated. We 
therefore performed fits on all 16 3 x 32 data (five 
k values) using a large variety of fit intervals (de- 
manding a common slope for the corresponding 
SS and SL correlators). We looked for fit inter- 
vals of at least three units of time which had x 2 
per degree of freedom near one. We selected the 
lowest value of the lower limit of the fit range that 
met these conditions. For the 16 3 x 32 data, the fit 
range so selected was 5-9. For the single k value 
at 12 3 x 24, we selected 3-9, and for 24 3 x 48, 
we selected 4-12. Intra-kappa correlations were 
included|8|. We will discuss the effect of varia- 
tions of the fit range and inter-kappa correlations 
in Subsection 4.3. 

4.2. Results for fs 

In Figures [| and |] the horizontal line superim- 
posed on the local effective masses is the common 
mass obtained from a simultaneous jackknife fit 
to the SS and SL correlators using the fit inter- 
val 5-9. The k=0.158 results are representative 



of the r 


esults for each of the five 


k values 


. Once 


Table 4 












Fitted parameters as a 


function of k at j3 


= 5.9 


K 


0.159 


0.158 


0.157 


0.156 


0.154 


Sb 


324 


316 


356 


344 


364 




11 


10 


11 


10 


10 


ma 


0.659 


0.662 


0.687 


0.687 


0.710 




0.014 


0.010 


0.009 


0.008 


0.007 


m'a 


0.940 


0.935 


0.950 


0.945 


0.954 




0.015 


0.013 


0.013 


0.012 


0.012 



the common mass and two intercepts have been 
obtained, the value of fs at each k value is deter- 
mined from the intercepts of the SS and SL fits. 
These results are presented in Table |[ which also 
contains the effective mass of the meson in lattice 
units as a function of n. Below each quantity is 
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its jackknife uncertainty, m! will be discussed in 
Section 5. 

We then extrapolate Jb to k c the critical value 
of k. We measured k c = 0.1597 ± 0.0001. Since 
the dependence on k is weak, the uncertainty in 
k c has a negligible (of order 1 MeV) effect on 
the results. The results for the extrapolation as 



380 
360 
340 
320 
300 
280 




154 155 156 157 158 159 160 
Figure 5. fs as a function of k. 



well as the results in Table [| are plotted in Fig- 
ure |. At k c , we find for the decay constant 
f B = 319 ± 11 MeV (jackknife). The slope with 
respect to kT 1 is 188 ± 32 MeV. The extrapolated 
value of fs and slope with respect to k _1 are su- 
perimposed on the Monte Carlo data. The value 
plotted at k=0.1597 is the extrapolated value of 
fs with its errors. The nonlinearity from using k 
as the horizontal axis rather than k _1 is imper- 
ceptible in the superimposed fit, and were we to 
have fit linearly in k rather than k _1 we would 
have changed the extrapolated value of fs and 
the (transformed) slope negligibly relative to the 
statistical errors. We note that the x 2 was 51 for 
44 degrees of freedom. 

We conclude this section by noting our results 
for fs in the static limit from (3 = 5.7 k = 0.165 
and (3 = 6.1 k = 0.151. In the former case us- 
ing the fit range 3-9, we find /b=351 ± 13 MeV 
(a^ 1 — 1.15 GeV) and in the latter case using the 
fit range 4-12, we find / s =359 ± 7 MeV (o -1 = 
2.43 GcV). Using the pion mass to determine the 
corresponding value of K for (3 = 5.9 we obtain 



from our fit to the f3 — 5.9 data the correspond- 
ing values of fs- They are Jb = 330 MeV for 
13 = 5.7 and f B = 367 MeV for (3 = 6.1. 

4.3. Dependence on Analysis Procedure 

In this subsection, we investigate the depen- 
dence of our results on (1) the fit interval, and 
(2) the inclusion of inter-kappa entries in the co- 
variance matrix. Perhaps the most interesting of 
the variations in analysis is changing the admix- 
ture of the excited state. This will be discussed 
in the Section 5. 

The dependence of Jb on fit interval has been 
investigated by rerunning the analysis over the 
other fit ranges, in addition to the primary 
range 5-9. We have reproduced the analog of 
Figure || in Figure |^. The results from the addi- 
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Figure 6. Jb function of k. 



tional fit ranges are displaced slightly from their 
true k values. From left to right, they are 6-10, 
5-9, 4-8, and 3-7. In general, the values are con- 
sistent with those of the primary fit range, 5-9, 
except at small values of k where it has already 
been noted that the \ 2 P er degree of freedom in- 
dicates that excited states are contributing when 
the lower limit of the fit range is less than 5. 

We now consider the impact of inter-kappa cor- 
relations. It is apparent from comparing the 
k = 0.158 effective mass plots with those of other 
k values (not depicted) , that there are strong cor- 
relations between corresponding quantities at dif- 
ferent n values. The reason these were not in- 
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eluded in our analysis is that with 50 data points 
per configuration (5 k's and 5 T's for SS and SL) 
it is impossible to compute a nonsingular covari- 
ance matrix. To investigate the effect on the ex- 
trapolated value of /s, we therefore reduced the 
number of k values to 2, selecting k=0.158 and 
k=0.156. We then performed the fit with and 
without the inter-kappa entries of the covariance 
matrix, and found that a decrease in the fitted 
and extrapolated values of fs resulted from in- 
cluding inter-kappa correlations. However these 
effects are in all cases less than the one sigma 
level. Larger statistics would allow us to study 
the impact of inter-kappa correlation with more 
k values in the simultaneous fit. 

5. EXCITED STATES 

In this section, we will examine evidence for 
the first radially excited state of the B meson, 
and study the impact of the admixture of this 
wave function. Preliminary results for the or- 
bitally excited states are presented in the second 
subsection. Finally, finite volume systematics are 
discussed in the last subsection. 

5.1. Radial Excitations 

We begin this subsection by looking at the ef- 
fect on the fitted values of fs of the admixture 
of the states created by the first radially excited 
state smearing function. The smearing function 
for the analyses in the preceding section was the 
linear combination of smearing functions which 
diagonalized a two-by-two matrix of SS correla- 
tors which was averaged over the same range as 
the primary fit range. For the 16 3 x 32 data, the 
admixture of the trial first excited state increased 
from 5% to 7% with increasing k. 

To see the impact of this 5 to 7% admixture 
on the final results, in Figure ^ we have replotted 
the data shown in Figure || along with the values 
of /b obtained from the undiagonalized smear- 
ing function. The additional data is displaced 
slightly from its correct k value for visibility. The 
extrapolated value of Jb using the undiagonal- 
ized smearing function is 333 ±12 MeV, a 4% 
increase. 

We now examine the linear combination of the 
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Figure 7. fs as a function of k. 



trial ground and first excited states which is or- 
thogonal to that used for the ground state. Since 
we do not have correlators involving states with 
higher radial excitations, we expect that this 
orthogonal state is missing contributions from 
higher radial excitations of the same order as the 
mixing of the trial ground and first excited state. 
It is neverthelesss interesting to study the effec- 
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Figure 8. k — 0.158 SS excited state effective 
mass. 



tive mass plots obtained from this orthogonal lin- 
ear combination, which is mostly the trial first 
excited state. In Figure || we plot the SS local 
effective mass for this correlator at K = 0.158, 
and in Figure || we plot the corresponding SL 
correlator. The same criteria to determine the 
fit interval as in the ground state were used, and 
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Figure 9. k — 0.158 SL excited state effective 
mass. 



gave a fit interval from 3-5. The excited state 
effective masses in lattice units, m'a, using this 
fit range are presented in Table 4. As usual, this 
is a common effective mass fitted to both the 55 
and SL correlator. The 25-15 splitting in lattice 
units, m'a— ma, decreases from 0.281 to 0.244 as 
k goes from 0.159 to 0.154. 

The decay constants of the excited state for the 
five k values are approximately 600 MeV with a 
variation in k and jackknife errors both less than 
±15 MeV. A slight trend is toward larger decay 
constants with increasing n. The value of the 
decay constant, f' B , extrapolated to k c is 618 ± 
11 MeV (stat). The quality of the plateau is much 
less convincing than for the (3 = 5.9 ground state. 
The two-state approximation may induce larger 
systematic errors here than for the ground state. 
In addition, finite volume effects may be larger 
for the excited states. One indication supporting 
these preliminary results is that the values of the 
effective mass and decay constant are increased 
by about one sigma only if the fit range is changed 
to 4-6. 

5.2. Orbital Excitations 

The lowest P wave heavy-light mesons have 
light quark total angular momenta ji = 1/2 
and 3/2. Each of these states is degenerate with 
the J = 0, 1 for the ji = 1/2 state and J = 1, 2 
for the ji = 3/2 state. Smearing functions for 
these P wave heavy-light mesons can be gener- 
ated from the SRQM with the mass parameter 



\x which gives the best fit for the S wave ground 
state. 

To date only a preliminary analysis has been 
performed on the P waves. The results for the 
mass splittings are given in Table 5. Considerable 



Table 5 

The mass splittings (in MeV) for the lowest 
heavy-light P states at j3 = 5.9 for various light 
quark k values. 



K 


IP1/2 - 15 


lP 3/2 - LP 1/2 


0.158 


386 ± 90 


43 ± 102 


0.156 


395 ± 86 


59 ±95 


0.154 


416 ±81 


62 ±90 



refinement will be required before the splitting 
between the ji = 1/2 and 3/2 states can be ob- 
served clearly. 

5.3. Finite Volume Corrections 

The systematic effects of finite volume are un- 
der study. It would be expected that these ef- 
fects are more pronounced for the excited states 
than for the ground state because the RMS ra- 
dius for excited states is larger than for the 
ground state. The excellent agreement between 
the SRQM wavefunctions and the measured wave- 
functions of the 15, IP, and 25 heavy-light states 
reported by Thacker Q| . allows us to estimate the 
effects of finite volume with our periodic bound- 
ary conditions using SRQM results. We calcu- 
lated the static energies using 48 Coulomb gauge 
fixed configurations at /3 = 5.9 for lattices of spa- 
tial size 12 3 , 16 3 , and 20 3 . Using the mass param- 
eter \i needed in the SRQM (obtained from our 
study of the heavy-light mesons on the 16 3 x 32 
lattices) we estimated the effects on various mass 
and wavefunction parameters at the other two 
volumes. The results are shown in Table 6. The 
typical variation between results at 12 3 and 16 3 
are 10% or more, while the variation between 16 3 
and 20 3 have dropped to a few percent. The 
validity of these model results is presently be- 
ing checked by a complete analysis of heavy-light 
mesons on the 12 3 x 24 and 20 3 x 40 lattices. 
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Table 6 

Finite volume effects in the SRQM. Variation of 
eigenvalues and properties of the eigenstates for 
the ground and radially excited states. Results 
are at f3 = 5.9. e„ denotes the eigenvalue for the 
lS(n = 0) and 2S(n = 1) states. 

Measure TF W 20 3 





0.900 


0.912 


0.911 


ei - eo 


0.251 


0.278 


0.277 


1*0(0)1 


0.193 


0.209 


0.206 


l*i(0)| 
< r 2 >V2 


0.239 


0.236 


0.236 


3.320 


3.01 


3.03 


< r\ >V2 


5.09 


5.50 


5.53 



6. CONCLUSIONS 

By using good approximate wave functions in 
a multistate smearing calculation, it is possible to 
control the systematic errors associated with ex- 
tracting the decay constant from nonasymptotic 
time. Our results are very encouraging, not only 
for the present calculations, but for the determi- 
nation of other B-meson decay parameters and 
form factors. All such calculations require the 
external meson to be in a pure eigenstate (i.e. 
on shell). By using the optimized smearing func- 
tions discussed here, this requirement can be met 
with a minimum separation between the B-meson 
source and the electroweak vertex, greatly im- 
proving the precision of the results. 

We find that at (3 = 5.9 and k — ^critical 

/B .319 ± ll ( stat)x^x( r ^)MeV(4) 

The efficacy of our smearing method is even 
greater at /? = 6.1 as can be seen by comparing 
Figs. 1 and 2 with Figs. 3 and 4 with the same 
physical volume. 

In addition to ground state properties, this 
method allows extraction of the properties of the 
lowest radially and orbitally excited states. Our 
results for these quantities are preliminary, and 
we expect more precise estimates to be obtained 
by further analysis. 

Further numerical studies are required to check 
the finite volume effects, determine the effect of 
using an improved action for the light quarks, and 



include additional light quark mass values at (3 = 
6.1 which will allow a better determination of the 
scaling behaviour of these properties of low-lying 
heavy-light states. 
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